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Abstract 

We propose a semiclassical approach to calculate multiparticle cross sections 
in scalar theories, which have been strongly argued to have the exponential form 
cxp (A _1 F(An,e)) in the regime A — > 0, An, e = fixed, where A is the scalar coupling, 
n is the number of produced particles, and e is the kinetic energy per final particle. 
The formalism is based on singular solutions to the field equation, which satisfy 
certain boundary and extremizing conditions. At low multiplicities and small ki- 
netic energies per final particle we reproduce in the framework of this formalism 
the main perturbative results. We also obtain a lower bound on the tree-level cross 
section in the ultra-relativistic regime. 
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1 Introduction 



Recently, the problem of multiparticle production in weakly coupled scalar field theories 
has received close attention. This problem has been initiated by the qualitative observa- 
tion [[[], ^ that in the |0 4 theory at tree level the probability of processes producing large 
number of bosons exhibit factorial dependence on the multiplicity of the final state. This 
dependence originates from the large number of tree graphs contributing to the process 
of multiparticle production: at multiplicity n ^> 1 the number of graphs is of order n\. 
At n ~ 1/A this factor is sufficient to compensate the suppression due to the smallness of 
the coupling constant, and the tree-level multiparticle cross sections become large. Much 
efforts have been made to understand how this behavior is changed by loop corrections, 
though a conclusive result on this issue is still lacking right now. 

So far, quantitative calculations have been performed mostly at, or near, the multi- 
particle threshold, where some perturbative techniques have been developed and exten- 
sively explored. The tree amplitude of transition from one initial virtual particle to n 
real bosons at rest (the 1 — > n process) can be computed either by summing Feynman 
diagrams or by using some appropriate classical solution |3|, f|, || and the result reads 

n-l 

< ee (0) = n\ (^\ 2 (1) 

(the boson mass in this formula and further is set to 1). The same methods have been 
applied for calculating the amplitude beyond the threshold or beyond the tree level ||. 
In the first case it has been found that when the final particles are non-relativistic, the 
tree amplitude is an exponent of the total kinetic energy of final particles, 

A^ ee {e) = A^ ee {0)e-l ne (2) 

where A^ ee (0) is given by eq.(|l|) and e is the kinetic energy per particle in the final state. 
The exponential fall of the tree amplitude beyond threshold in eq.(^) is not sufficient, 
however, to make the cross section small. The second result concerns loop corrections at 
exact threshold and reads that the leading-n contributions from each loop level (namely, 
the An 2 contribution from the first loop, X 2 n 4 from the second and X k n 2k from the A;-th) 
sum up to an exponent, so at not very large n, (when subleading on n contributions can 
be neglected, presumably at An 1) the 1 — > n amplitude at threshold has the form, 

A n (0) = < ee (0)e BA ™ 2 (3) 
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where B is a constant that depends on the number of spatial dimensions, 

dk 9 



In particular, in (3+1) dimensions (d = 3), the numeric value of B is 

B = -^(ln(7 + 4v/3)-nr) 

The physically interesting quantity is however not the amplitude, but the cross sec- 
tion, or transition rate. For the 1 — ► n process near threshold this quantity is easy to 
evaluate at n -C 1/A, having on hand the two results above. In fact, if e is small enough 
for the amplitude to be constant in the whole phase volume, the cross section is equal to 

a(E,n) = \A n (e)\ 2 V n 

where A n (e) = A^ ee exp(— ~ne + BXn 2 ), and V n is the bosonic phase volume. An impor- 
tant point to note is that at large n, V n has the exponential form except from the factor 
of l/n\, 

V n oc -^7 exp ( ^ fln-^- + l) + ^— -ne + 0(ne 2 ) ) (4) 



n\ \ 2 V ird 
and now it is easy to verify that the cross section is exponential, 

a(E,n) oc exp ( nln— -n+ — (ln— 1 + l] + (- — - - - ) ne + 2ReB\ 2 n ) (5) 
\16 2 V 7ra / \ 4 3/ J 

Though eq.(|) is valid only at small e and An, the form of a(E, n) strongly supports the 
hypothesis that in the most interesting regime An ~ 1, e ~ 1 the cross section is also 
exponential 0. 

a(E } n) oc exp (— F(An, e)) , e= (6) 

VA / n 

Moreover, there are indications j7| that the exponent F(Xn, e) is independent of the few- 
particle initial state (i.e. the cross section of 2 — >• n, 3 — >• n, etc. processes coincide, with 
exponential accuracy, with that of 1 —>■ n). The function F(Xn, e) is unknown, but some 
terms of its expansion at small An and e can be found from (|5]), 

An , . d A e _\ „ . ( d - 2 5\ , ^^^,22 



F(An, e) = An In An + - In — + 1 An + Ane + 2ReB\ 2 n z + 

16 2 V to / \4 3/ 

+ 0(A 3 n 3 ) + 0(A 2 n 2 e) + 0(Ane 2 ) (7) 
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The situation that has emerged here shows a complete analogy to the instanton-like 
processes at high energies. As in the latter case, the exponential form of the multipar- 
ticle cross sections is a strong argument in favor to the semiclassical calculability of the 
exponent F, but says nothing about the nature of possible calculation schemes. 

In this paper we propose a semiclassical method to calculate the multiparticle cross 
section. By using the coherent state formalism, we reduce the calculation to the problem 
of solving classical field equation with certain boundary conditions in the asymptotic 
regions t — ► ±00. The technique, as well as the boundary value problem are very similar 
to the those in the case of instanton transitions. The most essential difference from the 
latter case is that the boundary value problem for multiparticle production possesses 
only singular solutions. In particular, the field configuration defining the cross section is 
the one that is singular at one point t = x = and regular elsewhere in the Minkowskian 
space-time. Note that some other approaches utilizing singular solutions have been 
proposed recently for both multiparticle processes in scalar theories and instanton-like 
transitions |8], [| [TU], [TIL most being inspired by the Landau procedure for calculating 
the semiclassical matrix elements \12\. We emphasize, however, that in the framework 
of our formalism the singular field configuration is determined in a unique way by the 
boundary conditions and the structure of its singularity in the Minkowskian space-time. 

The field configuration with the required properties can also be found in a different 
setting. Namely, if one makes analytical continuation to the complex times and looks 
for solutions to the boundary value problem which is singular on some surface (in the 
simplest case the surface lies in the Euclidean space-time) and extremizes the transition 
rate over all possible forms of this surface , one obtains the same field configuration as one 
would find by solving the boundary value problem. Actually, this formalism sometimes 
appears to be simpler and will be applied for making quantitative calculations in this 
paper, which include reproducing eq. (HD at small An and e and finding a lower bound on 
the tree-level cross section in the ultra-relativistic limit of the final state. 

This paper is organized as follows. In Sect. 2 we derive the classical problem for 
the multiparticle cross sections. Sect. 3 is devoted to the tree amplitude at threshold. 
A simplified, entirely Euclidean version of the classical problem is presented and the 
perturbative result for the energy dependence of the tree amplitude near threshold is 
reproduced. In the opposite, ultra-relativistic, limit, we obtain a lower bound on the 
tree cross section. In Sect. 4 we consider the amplitude (with loops) at exact threshold, 
where the procedure for calculating the amplitude is derived from the general formalism. 
We reproduce the exponentiation factor coming from leading-n loop in the limit An <C 1. 
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Finally, Sect. 5 contains concluding remarks. 



2 General formalism 

2.1 The boundary value problem 

We consider the scalar field theory without symmetry breaking in (d + l)-dimensional 
space-time, 

s = jd d+ \ Q(^) 2 -^ 2 -^ 4 ) 

The quantity in interest is the total transition rate from an initial few-particle state to 
all possible final states having given energy E and multiplicity n. One writes, 

a(E,n) = Y,\(f\P E PnSA\0)\ 2 (8) 

/ 

where A is the operator that creates the initial state from the vacuum, S is the S— matrix, 
Pe and P n are the projection operators to states with energy E and number of particles n, 
and the sum runs over all final states |/). Different choices of the operator A corresponds 
to different initial states: for example, A = corresponds to the 1 — > n process (for the 
operator A describing the 2 — > n process see J7j). We recall perturbative calculations in 
J7| indicating that F does not depend on the particular choice of A, providing the latter 
is independent of A parametrically. Making use of this fact, we will calculate (|8|) for 
the operator A most convenient for our purpose. Namely, we choose A in the following 
exponential form, 

A = e^°) 

where j is some arbitrary number. 



Following the technique of (T^|, we derive the classical boundary value problem for 



a(E,n). Using the coherent state formalism |TJ|], one rewrites eq.([8]) in the following 
integral form, 

<t(E, n)= J db* k db k d£ dr] V<\> Vtf exp (- J dk blb k e iuJki+lv + iE£ + inr]+ 

+ 5,(0,0,) + B f (bl(j> f ) +B*(0,^ + B)M' f ) +iS[<f>] -iSM +J0(O) +J0'(O)) (9) 
In eq.(^J), -B's stay for the boundary terms, 
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Bf(K,4>f) = ~Jdkb* k b\e 2 ^ T f + JdkV2^btM^ kTf ~ \ j dku^ f (k)<P f {-k) 

where Tf is some final time moment (the limit Tf — > +00 is assumed), while <fii(k) and 
0/(k) are the Fourier transformations of the field in the initial and final asymptotic 
regions. 

It is easy to see that if one takes j to be of order 1 / V^A, then the integral (^) has 
the semiclassical nature in the limit A — > 0, XE, Xn = fixed, and one can expect that 
it is saturated by a saddle point, where 0, 0', b^ and are of order l/y/X, and £ and 
f] are of order 1. In this case, the cross section certainly has the exponential form (|6|). 
However, the assumption j ~ 1/ \/~X contradicts the requirement that j does not depend 
on A, which forces us to take j parametrically smaller, j ~ 1. This difficulty is obviously 
the consequence of the non-semiclassical nature of the initial state that contains few 
energetic particles, in contrast with the final state consisting of many soft ones. A method 



to overcome this difficulty has been suggested for the instanton transitions |L5[ [L6| and 
can be applied for our problem in an analogous way. In this method one evaluates the 
integral @ for j = a/VX, where a is some constant, in saddle-point approximation, 
and find the cross section in the exponential form a ~ e w , where W is of order 1/A and 
depends on a. After that one takes the limit a — *■ 0. The claim is that in this limit one 
reproduces the value of W at j ~ 1. 

Apparently, the weak point in this way of reasoning is the limit a —>■ 0: it is neither 
obvious that this limit is smooth, nor that it reproduces the amplitude with one initial 
particle. Moreover, the exponent of the cross section of the 1 — > n process contains 
contribution from loops (the term 2KeBXn 2 in eq. ©) that seems to be of quantum, rather 
than classical, nature. Nevertheless, it turns out that this nontrivial exponentiated loops 
are reproduced by the semiclassical calculations (see Sect. 4 below), which is a sound 
argument in favor to the hypothesis about the a — > limit. 

Therefore, we assume that j ~ 1/ \f~X and look for the saddle point of the exponent 
of the integrand in (|). The saddle-point equations can be divided into two groups. The 
first group contains the equations for (f>, 

g = ij5(x) (10) 

i0 i (k)+w k i (k) = O (11) 
i<j) f (k) - uj U (j) f (k) + V2^b\e l " kT ? =0 (12) 
- b k e iuJk ^ +iri - b\e 2luJkT f + V2^(f) f (k)e iuJkT f = (13) 



6 



and similar equations for 0'. The equations from the second group relate the energy E 
and the number of final particles n to other parameters, 



E= dku^lh^^ (14) 

n = Jdkblb k e iuJ ^ +iri (15) 
First, let us consider the equations for 0. Eq. ([!(]) is simply the field equation with a 



(5 like source, while eqs.(ll), (12), fllg) can be considered as boundary conditions in initial 
and final asymptotics t —>■ =poo. It is convenient to rewrite these boundary conditions in 
a more transparent form. For this end we note that since is a superposition of plane 
waves in the limit t — > — oo and t — ► +oo, eq.(|ll|) can be satisfied only if the initial 
asymptotics of is purely Feynman, 

<f>i{k) = a\e iuJkt , t -> -oo (16) 

where a-^ are arbitrary Fourier components, while eq.(|13D implies the following asymp- 
totics of in the final asymptotic region, 

; (k) = -== (6 k e^+ tr ™' + b\e luJkt ) , t -> +oo (17) 



V2^k 

It is easy to see that eq.(|T7D satisfies the condition flT2] ) automatically. 

Let us turn to the equations from the second group. The physical meaning of eqs.(|i~4T) 
and (0) is simple: they read that E and n are the energy and the number of particles 
of the field in its final asymptotics (fl7|) . Since satisfies the sourceless field equation 
at all values of t but t = (where the source is not vanishing), the energy of the field 
conserves in the two regions t > and t < separately, but may have discontinuity at 
t = 0. The energy in the region t > is equal to that in the limit t — >• +oo and therefore 
is £?. To find the energy at t < we make use of the 0's initial asymptotics, and since 
the latter contains only Feynman components, it vanishes. So, we find that the energy 
has a finite jump at t — which, naturally, is associated with the 5-functional source 
located at t — x = 0. 

To simplify further discussions, let us make two conjectures that, as we will see in 
what follows, do not lead to contradiction. The first is that the saddle-point values of 6k 
and 6 k are complex conjugated to each other. The physical meaning of this assumption 
is that the sum over final states in eq.(@) is saturated by a single coherent state. The 
second conjecture is that the saddle point values of £ and rj are purely imaginary, 

f = ~iT, r) = i9 
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where T and 9 are real (for further convenience we choose different sign conventions for 
T and 9). 

With the two conjectures formulated above, the field configuration describing the 
multiparticle process at given energy E and multiplicity n is the solution to the field 
equation with source 

g = v5 d+ \x) (18) 

with the two boundary conditions, 

0(k) = a k e^ fc *, t -> -oo (19) 

0(k) = - 1 = (b^T- 9 -^ 1 + b\e iuJkt ) , t -> +oo (20) 

It is easy to notice that the initial boundary conditions (|I9| ) can be reformulated in 
Euclidean language. In fact, making the Wick rotation to the Euclidean time r = —it, 
eq.(|T9"D reads that 0(k), as a function of r, contains only the decaying component in the 
asymptotics r — > +oo, 

0(k) = a k e"^ T , r -»• +oo (21) 
In contrast, the the final asymptotics (|20| ) contains both frequencies and cannot be rewrit- 



ten in Euclidean language. Moreover the field in the final asymptotics is, in general, 
complex. 

The boundary value problem for <p' can be derived in analogous way, 

% = 

(f)'(k) = a]l e- iu)k \ t -> -oo 

(j>'Ck) = —L= ( b k e~ tuJkt + b* e ^ T - e+Wkt ) , t -> +oo (22) 
V2^k v ' 

Notice that if is a solution to the boundary value problem (|P8"| , [T^ |20[), its complex 
conjugate 0* satisfies eqs.(p2|). This fact simplifies our calculations, since we need to 
solve only one boundary value problem instead of two. In further discussions we will 
assume that </>' = <fi* . 

The saddle point of the integral @ determines the cross section, which has the 
exponential form, 

a(E,n)~e w ^ 

where 

W{E, n) = \E{Xn, e) = ET-n9- 21mS[<f>] (23) 
A 



8 



The relation between E, n and T, 9, can be found either from eqs.(l4[) and (|15D, or, 
equivalently, from 

2 d_hnS = E 2 dj™s =n m 
dr ' dd n - [ ] 

which can be easily understood if one recalls that £ = —%T and rj = id are the saddle 
point of the integrand in eq.@. From ( p3|) one sees that W(E,n) is the Legendre 
transformation of 2ImS'(T, 9), and therefore one obtains the following important relations, 

dW dW 

Having derived the boundary value problem for calculating the transition rate at finite 
j, let us discuss the limit j — > 0. It can be shown that in this limit the field configuration 
becomes singular at x — 0. In fact, according to eq. (pT8|) , has discontinuity at t — 0, 

50(x) = 0(x)| t=+o - 0(x)| t= _ o = ij6 d (x) 

This discontinuity leads a jump of the energy at t = 0, since the latter contains the term 
| /o?x0 2 . On the other hand, the discontinuity of the energy E is supposed to be finite 
while that of <fi is proportional to j and tends to 0. When j is small one has 



E = Jdx 0(0, x)50(O, x) = ij<j>(0) 



One sees that when E is fixed and j — > 0, 0(0) goes to infinity, which means that the 
field configuration becomes singular at x = in the limit of vanishing source. This is not 
surprising, since in the limit j — > eq. (|I0D becomes the sourceless field equation, whose 
regular solutions conserve the energy and therefore do not obey the boundary conditions. 

So, to evaluate the transition rate one should find the solution to the field equation 
which obeys the boundary conditions and has singularity at t — x = 0, but remains 
regular elsewhere in the Minkowskian space-time. For doing calculations, however, we 
will use another formulation of the boundary value problem. 



2.2 Extremization procedure 

Let us discuss the structure of singularities of our solution. Recall that in Minkowskian 
space-time, is regular everywhere except a point-like singularity at x = 0. However, if 
one extrapolates into Euclidean times, it may occur that develops more singularities 
beside that at x = 0. Let us consider a simple possibility that is singular on some 
(i-dimensional surface A in the Euclidean space-time, which we will parametrize either 
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as x M = x^Si), where Si are d coordinates on the surface, or r = tq(x). In the region 
near A, is inverse to the distance to the surface, 




where l(x) is the distance form the point x to A. In fact eq.(^) is the only possible form 
of the leading singularity of in the region near the singularity surface A. In fact, one 
can even develop a perturbation theory on I and find the correction to eq.(p6|), which 
is of order l(x). However, one will see soon that ambiguity begins at the order of 0(7 3 ) 
(on other words, the terms higher than 0(l 3 ) is not defined uniquely by the form of A), 
which reflect the fact that the solution is not defined uniquely by the surface where it is 
singular. If two solutions are singular on the same surface A, the difference between them 
goes to when one approaches the surface as Z 3 . Note that if A touches the plane r = 
only at one point x — 0, then in Minkowskian space-time has the required structure, 
i.e. is singular only at x — 0. 

Let us take an arbitrary surface A satisfying the latter requirement and determine 
a field configuration 0, which consists of two parts <f>% and 02, as follows. Both 0i and 
02 are supposed to be solutions to the field equation and singular on A, but each of 
them obeys one boundary condition from eqs. (|19D and ( p0|) . Namely, we require that 0i 
decreases in the Euclidean asymptotics, 

01 — > 0, T — > + OG 

while 02 obeys the boundary condition in the Minkowskian limit, 

02 (k) = (b k e u " T - e - 4u " t + 6i k e^') , t -> +oo 

One may imagine the field is defined on a particular contour on the complex time plane, 
which at each value of x goes along the Euclidean time axis from zoo to some fro(x) lying 
on the singularity surface and then goes back to and then along the Minkowskian time 
axis to oo (fig.0). The field 0i is defined on the first part of the contour, (ioo,fr ), while 
02 on the two final parts, (fro, 0) and (0, oo). Note that there is a region on the Euclidean 
time axis where both 0i and 02 are defined, namely (iro,0). 

Despite the fact that obeys the boundary conditions of the boundary value problem, 
it is not the solution to the latter, since 0i and 02 need not to be equal at t = 0. In 
other words, for a generic surface A, contains discontinuities on the whole plane t = 0, 
instead of being singular at one point t = x = 0. On the other hand, if one manages to 
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choose the surface A in such a way that 0i(O,x) = 2 (O,x) for any x 7^ 0, would be 
the solution to the boundary value problem. 

Let us define the Euclidean action of as the sum of the action of 0i and 02, S[<f)i] 
and 5 , [0 2 ], each calculated along the corresponding part of the contour, 

T0(x) OO 

3 E [<I>] = - J dTdx(^(dM 2 + V{<h))- J dTdx(-(dM 2 + V(^-iJdtdxL(fo) 

+00 r ( x ) 

Since 0^2 are singular at r = r (x), both 5[0i] and £[02] are infinite, but since the 
integration contour for 0i and 02 goes at different directions near the singularities, one 
can hope that their sum S[<p] is nevertheless finite. 

Now we will show that the "correct" singularity surface A determined by the condition 
that 0i = 02 at t = can be found by extremizing the real part of the Euclidean action 
Se[4>] with respect to all possible form of the surface A, with the requirement that the 
point t = x = lies on the latter. First let us regularize to avoid dealing with infinities 
in intermediate calculations. For this end we replace the condition that is singular on 
A by the condition that = 0o on the same surface, where 0o is some large, but finite 
number, which will eventually tends to infinity. So, we set 0i = 02 = 0o on A. Since 
0i and 02 are, in general, different in the region near the surface, we expect that the 
derivatives of 0i j2 are different on A. Let us denote 

d n (</>i-<h)=j(si) (27) 

where d n is the derivative along the direction normal to A. The configuration 0, thus, 
can be regarded as the solution to the field equation with a source that is distributed 
over the surface A, 

dS E 



j(x) = / dsi j(si)5(x - Xp(s)) 

(in this formula we assume that the appropriate metric factor has been included in ds^. 
Now, as an intermediate step, we show that once 5e[0] has been extremized with respect 
to A, the source j(x) is proportional to 5(x) (in other words, the source is located at the 
point x = rather than distributed over A). Let us take an arbitrary singularity surface 
A and deform it slightly to A'. This can be represented as shifting each point x^ on A 
by a small vector 5x M = n^dx, where is the unit vector perpendicular to A, so that 
x^ + 5n M lies on A'. To ensure that x = is always a singular point, we require that 
6x^=0 = (fig.0). 

The new surface A' thus corresponds to new configurations 4>[ 2 = 0i,2 + <50i,2- We 
will evaluate the variation of 5[0i] and S[4>2] separately. The variation of 5[0i] is due to 
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two factors: the first is the variation of the field, 0i — » 0i + 50i and the second is the 
change of the integration region. The first contribution can be reduced to a boundary 
integral, since 0i is a solution to the field equation, 

T0(x) 



- JdrdxS QOVO 2 + V(<I>S) =- Jdx [dpfad^fa-VfaW 



ds {d^i-n^Sfa 



If is small, the second contribution that is associated with the change of the integra- 
tion region can be also reduced to a boundary integral, where each point on the surface 
A is integrated with the weight 5x. Therefore, the full variation of S[<pi] is 



5S 



E\<Pl 



ds 



Let us make use of the fact that 0i = 0o on A, and <f>' x = 0o on A'. We have 

5<p{x^) = (j>'{x^ - <p(x^) = (p'iXfj,) - <p'(x^ + Sxfj) = -{d^ ■ n^)5x(s) 
Therefore, 



6S[fa 



ds 



(^01 



ds 



-(dM 2 -v( ( f> 1 



8x(s) 



(2f 



where we have made use of the fact that on A the derivatives of along directions tangent 
to A vanish (since is constant on A). 

Let us turn to the variation of 5*[0 2 ]- The computation is completely analogous to 
the case of if>[0i], with the exception that now there is a boundary term at t — +00. So 
we have, 



5S 



Em 



ds 



(<9 n 2 ) 2 - V{(f> 2 ) 



- i J rfx<9 o 2 50 2 



t=+oo 



(29) 



However one can show that the boundary term at t — +00 is purely imaginary. In fact, 
since 02 and 2 obey the boundary condition ( ^0|) with the same T and 8, we have at 
t — > 00, 



1 



5bi L e 



io^T—O—iojjst 



5b* e 



t 



-00 



Substituting this, as well as the asymptotics for 2 , to the boundary term at t = +00, 
we see that the latter term in the r.h.s. of eq.(|29|) is in fact purely imaginary 

i J dx<9 o 2 <50 2 = J j0- (b k 6bt - bt5b k ) e^ T - e 
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and since we are interested only in the real part of the Euclidean action, this term can 
be dropped it in further calculations. Taking the sum of (p8|) and (|29|) , one finds, 



Now when one takes O to be large, both <9 n 0i j2 are large but the difference between them 
is small. Making use of eq.(p7|) one writes 



Now we see that the requirement that £ReS'[0] = for all variations A obeying Sx\ x= q = 
can be satisfied only if j(s) is proportional to the delta function, j(x) = joS(x). So, when 
we extremize the real part of the action varying the surface A, the source j, which 

is at first distributed over A, gather to a localized delta-functional source j(x) = joS(x) 

Suppose that we have performed this extremization procedure. Since now <fii and 
<p2 are equal on A and there normal derivatives are also equal (except from the point 
x — 0), these fields coincide in the region where they are both determined, namely, at 
r (x) < r < 0. In particular, <pi = <j) 2 everywhere on the plane r = but r = x = 0. 

So far we have been dealing with the regularized field configurations. Let us now take 
the limit 0o ~~ > 00 • The coincidence of <p\ and 02 at t = remains in this limit, however 
now the point t = x = becomes singular. The strength of the delta-functional source, 
jo, goes to in order to keep the jump of energy finite. So, (j) obeys the sourceless field 
equation in Minkowskian space-time, and is singular only at x — 0, thus it is the solution 
to the boundary value problem. 

To summarize, we have shown that the solution to the boundary value problem can be 
found by extremizing the real part of the Euclidean action Se[4>] (or the imaginary part 
of the Minkowskian action) over all singularity surfaces A containing the point t = x = 0. 
Let us now apply this formalism to find the cross section in various limiting cases. 

3 Tree— level cross sections 
3.1 General consideration 

Consider the exponent for the cross section, F(Xn, e), at small An. Keeping in mind the 
formula (0), one expect that at e ~ 1, An <C 1, the function F has the following form, 




ds ((dnfa) 2 - (dnfc) 2 ) 6X(S) 



A 





F(Xn, e) = An In An - An + An/ (e) + 0(A 2 n 2 ) 



(30) 
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where /(e) is some function of e. In what follows we will neglect the terms 0(X 2 n 2 ) and 
higher. Since these terms come from loops, it is equivalent to considering the tree level. 
Another way to see this is to compute the cross section corresponding to Q3"0|), 

a(E, n) ~ exp ( -F(nX, ej) ~ n!A n e n/ ^ (31) 



.A 

We see that the cross section depends on the coupling constant A as A n , which is natural 
for tree diagrams whose number of vertices is n/2. Even at tree level, the cross section at 
arbitrary e has not been calculated (for lower bound see ref. ||18|| ). In the non-relativistic 
limit e < 1 the perturbative result ([F]) yields the following formula for /(e), 

/(£) = -lnl6 + -ln- + - + (— -- 3 )e + 0(e>) (32) 

In this section we make no attempt to compute /(e) at arbitrary value of e. Our main 
goal is to reproduce eq. (|32|) from the formalism of Sect. 2. We will also try to estimate 
the tree cross section from below in the limit e — ► oo. 

Before considering the small-e limit, let us first point out that the calculation pro- 
cedure can be considerably simplified if one restricts himself to the tree level. We will 
work in the "extremization" formalism, not in the framework of the original Minkowskian 
boundary value problem, so we take an arbitrary surface A and calculate Se[4>] (we will 
deal only with the Euclidean action, so for simplicity further we will drop the index E). 

First note that from eqs.(EF) and fl3T)|) one finds 



A dn 16 J W 

So, in the limit Xn —>■ 0, 6 ^> 1 independent of the form of the function /(e). We see that 
in the final asymptotics the Feynman part of 0, b ke ~' luJkt+UJkT ~ d is much smaller than the 
anti-Feynman part, b^e luJkt . 

Recall that the initial asymptotics of <\>\ is the same as that of eq.(^TJ). By con- 
struction, 0i is defined on the first part of the contour of fig.|l|, (ioc, 0), however one can 
always analytically continue <\>\ to other parts. Let us investigate the difference between 
0i and 02 on the parts 2 and 3. Denote 



so is regular (in fact, equal to 0) on the surface A and obeys the boundary condition 

1 



0(k) = —L= (b ke -^ kt+UJkT ~ e + (6* k - a\)e lLOkt 
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in the final asymptotics. Let us show that the anti-Feynman part of 0, 6!_ k — al k is 
small and proportional to e~ 9 . Supposing that this is true, then the whole is a small 
perturbation on <pi and thus obeys the linearized equation 

(aJ + l + 3A0?)0 = O (33) 

with two boundary conditions 

4>U = o 

4>(k) = —^= a^-^t+^T-e + &ny ant j_Feynman part, t -> +oo (34) 

(we have made use of the assumption that 6 k ~ a k) • Since the equation is linear and the 
final boundary conditions contains a factor of e~ e , the solution is also proportional to 
e~ , which is consistent with our starting assumption. 

Now making use of the smallness of 0, the action, up to the contributions of order 
e~ d , is equal to 



S[<f>] = Si[0!] + % 3} [0i] " / dx [(^00(^0) - V'{<k)4>\ (35) 

J {3} 

It is easy to see that the action of 0i on the whole contour (the first two term in the 
r.h.s. of eq.(^)) vanishes (one can see this, for example, by making the Wick rotation of 
the third part of the contour). The last term in eq.(|35|) is reduces to boundary integrals. 
The boundary term on A is equal to since near A 0(x) tends to as I 3 , where I is the 
distance from x to A, while <90i ~ l~ 2 . The boundary term at t — +oo is 

S[<f>]=i 1^x0^0x1^+00 = ^Jdkala^- 9 (36) 

So, to compute the action, one need not really solve eq.(^): only the knowledge of 0i 
(more precisely, its Euclidean asymptotics) is required. Note that the Euclidean action 
is real. 

Once the action is found, one should extremize ( p6[ ) over all solutions that are singular 
at x = 0. In fact, one can see that one should extremize only fd\tal_ai L e UJkT , since e~ e is 
just an overall factor in eq. (|3"6]) . Moreover, since the action is bounded from below by 0, 
the extremum of the action is most likely the true minimum. We also expect that the 
action contains a classical factor 1/A, so, the extremized action has the following form 

2S(T, 6) = i e ^- e (37) 
A 
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where g(T) is some function. Eqs. (|24|) then read 

E = 2 w = j^ T i e, ' T> ~" < 38 > 
„ = _ 2 | = (39) 

These equations should be solved with respect to T and 6. Dividing to fl5P|), one 
obtains 

l + e = </(T) (40) 

We see that the parameter T depends on e but not on 9. Regarding ( flOf) as an equation 
on T, we denote its solution as 

T = T(e) 

and from eq.(|3£f) one finds 6 as a function of e and n, 

9 = g(T(e))-H*n) 
Let us substitutes the solution to the exponent of the cross section. One obtains 

\\V = ET - n6 - 2S(T, 6) = n(l + e)T(e) - n{g{T(e)) - In Xn) - n 
A 



Comparing the last equation with (p0|), one finds 

f(e) = (l + e)T(e)-g(T(e)) (41) 

Therefore, the problem of finding the tree cross section at any value of E and n 
can be formulated entirely in the Euclidean space-time. One looks for all solutions 
0i(t, x) to the Euclidean field equations which are singular at r = x = and decay 
as t — > +oo, and calculates for each solution the corresponding Fourier components 
from its asymptotics at r — > oo. Then one should maximize the integral fd\t alxi]f L e u)kT 
and determines the function g(T) (eq. (p7|)). The required /(e) can be calculated using 
eq. fllTD , where T(e) is the solution to eq. (|40|) 

Unfortunately, we are unable to carry out this program for arbitrary values of e due to 
the non-linearity of the field equation. At small e (which corresponds to non-relativistic 
final states) we will see that it reproduces the result found by perturbative calculations. 
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3.2 Non— relativist ic regime 
3.2.1 Leading order 

At small e, the typical momentum of final particles is kg = e 1//2 of much smaller than the 
mass, so one can expect that <fi\ is a slowly varying function of x. The x-independent 
solution to the field equation which is singular at r = and decays as r — ► oo is known 
explicitly, 

*M = &L- (42) 

Let us restrict ourselves to the leading order. The field configuration in the first approx- 
imation can be obtained from (f|^) by a simple modification, 

0x (r, x) = J \ 1 (43) 

V A smh(r — t (x)) 

where t (x) is a slowly varying function of x. One can check that eq.(f43"D satisfies the field 
equation to the accuracy of 0((<9 x To) 2 ). Note that r = r (x) is the surface of singularities 
of 4>i. The point r = x = should lie on the latter, so we require that r (0) = 0. 

First, from eq . (|43"D one can relate to tq(x). At r — > oo, the asymptotics of <j>i in 
eq.(g3} is _ 

0(r,x) = W^e ro(x) e- r (44) 
V A 

On the other hand, <f>\ can be expanded into plane wave in this region, 



(27T) d / 2 v^ 

Recalling that typical momentum k is small, as the first approximation one can replace 
co>k by 1 and the r.h.s. of eq. (|45|) behaves like e~ T , as that of eq . (|44]) . Comparing eqs.(^) 
and P5|), one obtains, 

,/k 0k e ta = */^e'»W (46) 



{27l)d /2 — y A 

Taking in this equation x = 0, and recalling the requirement r(x) = 0, one finds the 
following constraint on a^, 

(2^ ak = VT (47) 

To find one should extremize the r.h.s of eq. fl3"fi"D with the constraint (^7|). We 
follow the standard technique and introduce the term 
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where C is the Lagrange multiplier to the r.h.s. of eq.(|36"D and, varying with respect to 
a k , we obtain a k , 

a k = Ce~" kT (48) 



It is easy to notice that fl4q) in fact the minimizes of the action. The constant C can be 
determined from eq.(f47|), 

C = ^T d / 2 e T 
V A 

Let us now calculate the function g(T). According to eq. (|37|) , 

g(T) — laxfj dk a k a k e WfcT ) = T+^ ln(2vrT) + In 16 
so eq.pOD becomes 

l + ^ = l + e 
IT 

whose solution is obviously 

T(e) = ± (49) 
Now substituting T(e) and g(T) into eq.(|4"T|) one finds finally 

/( C ) = ^-^ln^-lnl6 
JK ' 2 2 e 



which is nothing but the leading terms in eq.floz]). 

Before considering the 0(e) correction to /(e), let us discuss the form of the surface 
of singularities r = r (x). From eqs. fl46|) , (f48|) one finds 

r (x) = = |x 2 (50) 

One sees that the surface of singularity has the form of a paraboloid, whose curvature is 
proportional to the typical momentum of the final particles k . At small e the curvature 
radius is much larger than the inverse boson mass, i.e. 1. One should keep in mind, 
however, that eq. (|50|) is valid only for not very large values of x (not much larger than 
k^ 1 . At larger x the actual behavior of the singularity surface is unknown. 

3.2.2 Next— to— leading order 

To find the first correction to the solution, eq.([|3|), let us substitute the ansatz 



A sinh(r — t (x)) 
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to the field equation. We expect that 4>i is suppressed, compared with the leading order, 
by a factor of k$, where ko is the typical momentum of the final particle, so we drop 
terms containing <fi\ and <p\. We expect also that each derivative with respect to x adds 
an additional factor of k . Denoting y = r — r , the equation for <pi is linear, 



Ocf) = 



sinh 2 y 



2 coshy 2 



sinh 2 y X \ sinh y sinh 3 y J J 

(51) 

We have ignored, for example, <9 x 0i, since it is of order 0(&q). The simplest way to solve 
eq. (|5lD is to make use of the following simple relations, 

Ocoshy = —6 °^ 2 ^ , Osinhy 



O- 



sinh y sinh?/ 
1 4 



sinh y sinh y 
to see that one solution to (1511) is 



U^ coshy+ (Ml smhy+ W (52) 




A y 6 6 2 sinh y 

However, one can also add to (|5^) an arbitrary solution to the homogeneous equation 
0(f) = 0. The general solution to the homogeneous equation depends on two arbitrary 
functions of spatial coordinates, -Fi^x) and has the form, 

= F 1 (x)/ 1 (r,x) + F 2 (x)/ 2 (r,x) 

where 

cosh y 
/i(r,x) = . 2 
sinh y 

cosh ?/ sinh ?/ 

sinh 2 ?/ 3 sinh y 

Note that fi has a double pole at y — r — r(x) = 0, while / 2 is regular on the singularity 

surface. As we expect that the field has a pole of the first order at y = 0, so we should 

rule out fi and set F\ (x) = 0. The function F 2 (x) can be chosen from the requirement 

that 0i decreases in the Euclidean asymptotics r — > +oo. In fact, in this limit ( |52D 

becomes 



/ 2 (t,x) = y— — — — (53) 



/|4(^r„ + (W) 



e r-r 



1 Actually, the term proportional to f\ arises when one shifts the singularity surface tq — > tq + 5to 
and try to expand the function sinh _1 (r — To — Stq) on Stq. Since we want the singularity surface to be 
r = to(x), this term should be excluded. 
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while the asymptotics of /2 is also growing 



/ 2 ( r ,x)~--e^-° 



but they cancel each other when we choose i^Cx) 



l /2r^2 



2 V A 



[^x r o + (^xTo) 2 ]. Therefore, the 



solution to the field equation with singularities at r = t (x) and decays at r — > +00, up 



to the order of k%, is 



sinh?/ 



+ 



6 



cosh y + 



6 



sinh y + 



(<9xT ) 



2 sinh y 



<9 2 r + (<9 x t ) 2 / coshy 1 . 
+ - I V — 2 . ~ 77 smh l/ 



(54) 



2 ^sinh^y 3" a sinh?/, 

To find at, in complete analogy with the leading order, one calculates the asymptotics 
of (f>i from eq. ([>§), 



2 - ^ 2 r ) + («9 2 r + («9 x t ) 2 ) ye"* 



(55) 



Let us note that, in contrast with the leading order, there is a term proportional to ye~ y 
in the asymptotics of <f>\. This structure also emerges when one evaluate the r.h.s. of 
eq.(f4"5|) to the next-to-leading order, replacing e 1 ^ 1 " by e _r (l — ^r). Introducing the 
function 



ax 



dk 



(27T) d / 2 v^ 



ikx 



eq. (f45j) reads 



h(r, x 



dk { 



r+ ' kx - ' a(x) + ^<9 2 a(x) ) <>- ; (5(5) 



Comparing eqs.(^) and (|56|), one finds the relation between a(x) and t(x), the next-to- 
leading-order version of eq . (|46D , 



ax 



I (2 - ^ 2 r - r (<9 x T + («9 x t ) 2 ) ) e*>W 



(57) 



Let us find out the constraint on a(x) comes from the requirement that To(0) = 0. Taking 
x = and noticing that <9 x To| x =o = since the singularity surface is tangent to the plane 
r = at x = 0, it is easy to show that this constraint is 



«( x ) + yr^ 2 a(x)| x=0 
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In momentum representation this relation reads, 

dk ( 5 



The difference of the constraint from that of the leading order is there is a small correction 
T^k 2 . The calculation is now straightforward and similar to that of the leading order. 
Introducing the Lagrange multiplier and taking variation with respect to a^, one finds 

a k = —j= 1 k 



V 12 

where C is a constant that can be found from ([58]). The final result of the calculations is 

d /, e \ /d-2 5\ 

/M = 2("^ + 1 ) + ' nl6+ (— -3j' f 

which, as expected, coincides with the perturbative result for the tree cross section near 
threshold. One can, in principle proceed further in this direction and calculate more 
terms in the expansion of / on the small parameter e by the same technique. In this 
way one could find 0(e 2 ) and higher corrections to eq. fl3"2]) , which have not been found 
by standard perturbative techniques. However, let us stop here and turn to the opposite 
limit e > 1. 

3.3 Ultra— relativistic limit in (3+1) dimensions. Lower bound 
on tree cross section 

In the ultra-relativistic limit e — > oo, the boson mass can be neglected. We will consider 
the most interesting case d = 3, where A is dimensionless. In this case the theory 
becomes scale invariant and the exponent of the cross section should be independent of 
the energy (the dependence of the pre-exponent on the energy can be determined solely 
by dimensional analysis). In other words, we expect that the function F(Xn,e) becomes 
independent of e at large e, 

lim F(Xn,e) = F(Xn) 

e^oo 

At tree level, this statement means that /(e) has a limit when e — > oo. 

From the first identity in eq.(^) one sees that T should tends to faster than 1/e 
(otherwise W = JdET(E) diverges in the ultraviolet). Now let us take the limit e — > oo 
in eq.(fnp. One finds, 

/(oo) = -g(0) 
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Recall that X~ 1 e 9 ^ is the minimal value of /<ika k cik over all solutions of the massless 
theory that are singular at r = x = 0. We are unable to find g"(0). Instead, we will try 
to bound g(0) from above. We will take a solution that is known analytically, calculate 
for it at and plug the result into the definition of g(0). If the extremum over is the 
true minimum (recall our discussion on the natrue of the extremum in subsection 3.1), 
this procedure would gives us an upper bound on g(0). 

We take the following trial singular solution to the massless field equation 0, 



P 



A x 2 + (r — p) 2 — p 2 
The Fourier components of this configuration is 

V A y u k 

which implies the following bound on g(0), 

g(0) < ln(87T 2 ) 
The corresponding lower bound on the cross section is 



er„, > n! 




Note that this lower bound grows factorially as n — > oo, which reflects the fact that at 
tree level not only amplitudes, but also cross sections becomes large at n ~ A -1 . 



4 Loop corrections at threshold 

4.1 General consideration 

In this Section we consider another limiting case. Namely, we take arbitrary An but small 
e <C 1, and will be interested in the exponent of the cross section to the leading order of 
e. In this limit, we expect from eq.([7j) that the result has the following form, 

F(An, e) = An In - An + ^ fin — + l) An + g(Xn) (59) 
16 2 V ltd J 

where the first term in the expansion of g{Xn) at small An should be 2ReSA 2 n 2 . The 
meaning of the function g(Xn) becomes clear when one recalls that in the limit e — > 
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the cross section is the product of the square of the threshold amplitude and the phase 
volume. Dividing the cross section corresponding to eq.(|59|) to the phase volume at small 
e, one finds the absolute value of the threshold amplitude, 

/ , \ n/2 

|A»|ocn!l-J e^ {Xn ~> 

Comparing with eq.(|I|), one finds that g is nothing but the contributions of the loops to 
the amplitude at threshold. Note that the reproduction of of even the leading term in 
g, 2KeBX 2 n 2 , is a nontrivial argument in favor to the validity of the whole semiclassical 
ideology, since it shows that the exponentiated part of loop contributions is semiclassically 
calculable. 

To find the form of the singularity surface, let us begin by recalling the result of 
Sect. 3 that e —>■ corresponds to T — > oo (eq.([EJ)), which means that the surface of 
singularities has a very large curvature radius and in the limit e — > can be considered 
as a plane. However, the discussions in Sect. 3 is based on the assumption that the limit 
An — > is taken. When one drops this assumption, one could expect that the presence 
of the source at x = deforms the surface of singularities near its location. This change 
should be local and the curvature of the singularity surface should tend to at large x. 
At finite An, one expects that the surface of singularities r = r (x) has the form similar 
to that sketched in fig|| The requirement of zero curvature at large distances from x = 
can be satisfied if To(x) tends to some constant as x — > oo. 

If the singularity surface is just the plane r = t^, the solution would be x-independent 
and equal to 

1 1 

A sinh(r — r a 

In the case when the singularity surface has the form shown in fig.^], the general form of 
is 

^2 1 



x • , , ^ + 0(r,x) (60) 
A smh(r - r^) 

where 0(r, x) vanishes at x — ► oo and decays into plane waves in the asymptotics regions, 
since it comes from the local (about x = 0) deviation of the singularity surface from the 
plane r = r^. 

Let us consider the asymptotics of at t — > +oo. Strictly speaking, the x-independent 
part of remains non-linear at any t and does not decays into plane waves. This is an 
artifact of the limit T — ► oo: at any finite T the field is a linear superposition of waves 
at sufficiently large t. To make linear, let us take the final part of the contour to form 
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at a small angle with the Minkowskian time axis. Namely, we take t = (1 + i5)a, where 
5 is small and a is a real parameter going to infinity. This leads to the replacement 

-> 2e (T ~ Too) 



sinhfr - r^) 

in the asymptotic region. This part of is homogeneous on x, so its Fourier transforma- 
tion is a delta-function in the momentum space. In contrast, the Fourier components of <p 
in the asymptotics t — > +oo are supposed to be well behaved functions of k. Denoting / k 
and (?k to be the Feynman and anti-Feynman components of <fi at the final asymptotics, 



4>(t,k) 



one finds that the asymptotics of is 

1 



/ k e-^ + #_ k e^ 



0(t,k) 



/ke — + ( 27r ) 



d/2 



V2^k 

Comparing eqs.fl2~Cf) and (|B"TD, one finds, 



-5(k)e To ° + £_ k e 



(61) 



(62) 



(2tt 



-5(k)e roc + # k 



(63) 



From eq.(^) one obtains 6 k = / k e~ WfcT+6 '. Since T is large, 6 k looks like a delta-function 
in the k-space: 6 k is non-zero only in a small region near k = 0, namely, k ~ T -1 / 2 . As 
/k is a smooth function of k, one can replace / k by its value at k = and write 



foe 



-uj k T+e 



Consider now eq.flBlf). Since 6 k has the form of the delta-function, vanishes. This 
important result means that </>(£, x) is purely Feynman in the final asymptotics. Eq . (|63|) 
now reads 

/ 1 R 

f* e -^T+e = ( 2vr )^. ^-5(k)e T °° (64) 
V A 

This equation should be understood as a symbolic representation of a relation between 
fo, T, 9 and r^. To find the latter let us take the integral of both sides of eq . (|64]) over 
dk. We obtain 

r L6 



fo 



j.d/2 e T-e+T 0l 
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The multiplicity and the (small) kinetic energy of the final state can be also expressed 



in terms of T, 9 and r^, 



l = J dkblb k e^ T ~ d = ^(2nTf 2 e T - e+2 ^ (65) 
/ dk y fc^e^ = ^(2nT) d ^e T -^ ■ ± ( 66 ) 



ne 

Solving eqs.(|65D and (B6l) with respect to T and 9, one finds 




^ = -ln^ + ^[ln- + l|+2r 

Note that T and 9 depend on the surface of singularities through but are insensitive 
of the behavior of the function tq(x) at finite x. In contrast, the action depends on the 
the precise form of this surface. In Sect. 3 we have seen that a part of the action (that 
is linear on 0) can be reduced to the boundary term <fid t (fi\t=oo and is equal to |. It 
is convenient to separate this contribution from the action and introduce the notation 
S' = ImS — S. The exponent of the cross section is now equal to 

W = ET — n9 — 2lmS = y (ln^- + ?j + n\n^-n- 25"[r (x)] - (67) 

where S' a functional depending on tq(x). Notice that the first three terms in eq. (|67j) 
coincide with those in eq. (|59D , we conclude that the exponent for the loop corrections to 
the amplitude at threshold can be obtained by maximizing the expression 

jg(Xn) = -n^ - S"[r (x)] (68) 

with respect to all possible surfaces of singularities To(x). The problem has a simple 
geometric interpretation: it is equivalent to finding the equilibrium configuration of a 
surface under the force n acting to the point x = 0, when the energy of the surface 
depends on its form through the functional S"[r (x)]. 

Let us summarize our results. To find the amplitude at threshold, one should 

1. Find for any function r = t (x) (which goes to a constant value r M as x -> oo) the 
solution to the field equation which is singular on the surface r = To(x) and decays 
at t +oo and at t = +oo has the form ( |60D where is purely Feynman. 
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2. calculate the action of the field, and find S' = ImS — | 




3. Extremize the r.h.s. of eq. (|68|) over all surfaces r = ro(x). 



Unfortunately, 5" seems to be a very complicated functional of To(x). One region 
when it can be calculated is the regime tq <C 1, where one can develop perturbation 
theory on r . Let us see that this is the regime of low multiplicities, n<l/A. 

4.2 Low multiplicities, An <C 1 

When tq(x) is small the field configuration can be found by solving the field equation 
separately at r 1 and r > r^, and matching in the intermediate region r C r < 1. 
For convenience we will deform the contour as shown in fig.f|, where r c is some value in 
the intermediate region, T ro < r c < 1. So, the contour in the t plane consists of four 
parts: (I) (ioo,ir c ), (II) (ir c , 0), (III) (0,rr c ), (IV) (ir c , ir c + oo). The energy on parts (I) 
and (II) is 0, while on parts (III) and (IV) of the contour it is equal to E. 

First consider parts (II) and (III). In these part r < 1, and the field is given by the 
following formula, 



Eq.(|69"l) is valid on both parts (II) and (III). While the leading singular terms coincide, 
there is difference between <p on (II) and (III). Denote 



Since we expect that this difference is much smaller than the leading term in eq . (|69D , 5<f) 
satisfies the linearized equation 



Furthermore, at r < 1, the mass and the spatial derivatives can be neglected, so 5(f> 
satisfies the equation (d^ + 3X(f>l)5(f) = 0, the general solution to which that vanishes on 
the singularity surface is 




(69) 



5(p(r, x) = m (r,x 



.) - (f) n (t, x) 



(d 2 T + dl 



+ 1 + 3X(f> 2 )5(f) = 




(70) 



where the function H^(x) is yet to be determined. 
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To the linear order on 50, the action on the parts (II) and (III) can be reduced to the 
boundary term at r = r c , 

S'jj+m = ~ jdxd T< t> ■ 50| T=Tc = j ((r c - Tqo ) J dx W(x) - J dx W(x)c(x)) (71) 

where c(x) = r (x) — r^. The next correction to S' II+III is suppressed by a factor of 
compared to the leading result (|7l"D and will be neglected. 

Consider now the parts (I) and (IV). On these parts of the contours, the solutions 
to the field equation can be represented in the form of the perturbative series on the 
background of the homogeneous field, 



+ 0(i) + 0(2) + • ■ ■ (72) 
where 0( O ) is the x-independent solution to the field equation, 

(73) 



A sinh(r - r^) 

0(i) is the linear wave on the background fl73[) satisfying the equation 

^0(i) - H0 (O ))0(i) = (74) 
everywhere except the singularities and having the following explicit form, 

I Ij^y 1°^^ on part W (75) 

xlwy l^ftM***' on P art ( IV ) ( 76 ) 

where Ck are some function of k, and /i 2 ( r ) are mode functions on the background field 



(i) 



(i) 



f}(r) = e-(— ) (ul - ^-J^ - - + 2 + : , , 



tanh(r — r^) sinh (r — r^) / 

f k 2 {r) = f^2 Too - r) 

and 0(2), etc. are higher corrections. The small parameter governing the expansion ( [72] ) 
is in fact r^. 

The relation between c k and the form of the surface of singularities can be found from 



the matching condition between eqs.(|72|) and ( p9|) at intermediate values of r. One finds 



c k = y dxc(x)e ikx 
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One can also relate c k with the function W(x). In fact, at r = r c < 1 one finds from 
eqs.(^|) and fl7|) _ 

' 1(Tl ~ T ^ ] * 1 ' /k W k c k e ikx (77) 



(i) 
where 



A 45 J (2n) d 



^ = 2^(^-1)^-4). (78) 

On the other hand, since r c r eq.([70D can be expanded on the parameter r /r c as 
follows, 

H = ^|w(x)(r c - r^) 3 - 3^(x)c(x)(r c - Too ) 2 + • • • (79) 

The r.h.s. of eq.(|77|) should coincide with the first term of the r.h.s. of ([79]), so we obtain 
the following relation, 

Let us now evaluate the action on the part (I). Up to the second order of the small 
parameter the action has the form 

oo 

Sj = Sj[0 O ] + / (^(0)^0(1) + ^0(0)^0(2) + ^(d^ {1) ) 2 + ^ , (0 ( o))(0(l) + 0( 2 )) + 

V(0o)(0(i)) 2 ) (80) 

Making use of eq.([74]) and the fact that 0(o) satisfies the field equation, the integral in 
eq. (|80D can be taken in part and the result is 

Si = Si[(f> ) - Jdx ^(0(i) + 0(2))<9 T 0( O ) + -0(i)9 r 0(i)^ \ T=Tc 

The action on the part (IV) can be calculated in an analogous way. Since S/[0o] + 
S/v[0o] — 0, the sum of the action on parts (I) and (IV) is given by the following 
boundary terms, 

S' I+IV = [dx (d T (o) ■ 54> {1) + <9 T (O) ■ 5<p {2 ) + \d T <p {1) ■ 50 {1) + \d T 5<p {1) ■ (1) ) (81) 



The analytical expressions for 02 are rather complicated. Fortunately, the calculation of 
the action requires only 502 at r = r c , which is equal to the second term in the r.h.s. of 
eq©, 

50 2 (r c ,x) = -3W(x)c(x)(r c - r^) 2 
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Substitute this to eg. (|STD one obtains 

S' I+IV = - (-(r c - roo) J dx W(x) + lJdx W(x)c(x)) (82) 

The full action can be obtained by taking the sum of fl8"2| ) and (|7T|). The dependence on 
r c disappears, as one can anticipate, and the action is quadratic on Ck, 

S'=y m) c(x) = ±J * (83) 



Now to extremize (BSD we note that 



so 



^-^"/(•^ (84 » 

1 /" dk 1 /• dk , l9 . , 

A S = n 7(2^ Ck -9A/(&p^ |Ckl (85) 



Differentiating g with respect to Ck, one obtains 

9A r 
~2~Wi 



9A n , , 

c k = — — (86) 



and 

1 9A 2 /■ rfk f 2 

— o = — n I - — —j = B\n 

\ y 4 J (2vr) d W k 

which coincides with the perturbative result for exponentiated leading-n loops. 

Let us find the condition for the approximation we use here to be valid. The small 
parameter is t^, and from eqs.(|84]) and (jB6|) we see that r ro ~ An. So, our calculations 
are reliable when n is small, n 1/A. 

Another remark should be made on the nature of the extremum. Since does not 
have definite sign (it is negative at k 2 < 3 and positive at k 2 > 3, see eq. (|78|)), the 
extremum of the r.h.s. of (|85|) is neither maximum nor minimum, but rather a saddle 
point (in contrast with the maximum for tree-level cross section considered in sect. 3). In 



the theory with broken symmetry is a positively defined function pl[ , and we have 
in this case the true minimum. 

5 Conclusion 

We have seen that the problem of calculation multiparticle cross sections can be reduced 
to a certain problem of the classical theory. The field configuration describing these 
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processes is a singular solution to the field equation with certain boundary conditions 
that optimizes the transition rate. Though most results obtained in this paper can also 
be found by making use of various perturbative methods, our discussions show that they 
can be derived from a single approach. We have pointed out an important fact that 
the exponentiated loop corrections can be calculated semiclassically. We also obtain a 
new result, namely, the lower bound on tree cross section in the ultra-relativistic regime. 
Hopefully, the formalism developed in this paper can be used in further analytical or 
numerical calculations of the multiparticle cross section 
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Figure 1: The contour. 




Figure 2: The deformation of the singularity surface. 
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Figure 3: The surface of singularities for configurations describing processes at exact 
threshold. 
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Figure 4: The contour for calculating the action at low multiplicities. 
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